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Extraction of the Spin Glass Correlation Length 
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The peak of the spin glass relaxation rate, S(t) = d[ — Mtrm^, t w )/H]/d£nt, is directly related 
to the typical value of the free energy barrier which can be explored over experimental time scales. 
A change in magnetic field H generates an energy E z — N s Xf c H 2 by which the barrier heights are 
reduced, where Xfc is the field cooled susceptibility per spin, and N s is the number of correlated 
spins. The shift of the peak of S(t) gives E z , generating the correlation length, £(t,T), for Cu : 
Mn Qat.% and CdCr1.rIno.3S4,. Fits to power law dynamics, £ (t, T) tx t a ^ T ^ and activated dynamics 
£(t, T) oc (int) 1 ^ compare well with simulation fits, but possess too small a prefactor for activated 
dynamics. 

PACS numbers: 75.50.Lk, 75.10.Nr, 75.40.Gb 
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The study of the irreversible behavior of the spin glass 
magnetization under a change of magnetic field allows 
exploration of the available states of a random frustrated 
system. There are various representations for the 

long time evolution and the dynamics of spin glasses, 
|^-|| but a coherent, overall accepted real space descrip- 
tion remains lacking. || The purpose of this paper is 
to extract a time and temperature dependent spin glass 
correlation length from a specially structured set of ex- 
periments, and to compare our results with available the- 
oretical predictions. 

The definition of a correlation length for a spin glass is 
difficult to express in measurable terms. Marinari et al. 
and Kisker et al. || introduced the time dependent 
equal time correlation function at time t. In the notation 
ofRef. 7, 



G(x,t) = V 



< Oi+xTi+xViTi >t 



(1) 



where the average is done at time t, and <jj (oi +!r ) and 
i~i (Tj+x) represent the z component of Ising spins at 
sites i (i + x) in two thermalized configurations in a 
box of volume V . To avoid accidental contributions to 
G(x,t), the two configurations are chosen to have zero 
overlap, q = V^ 1 J2i (J i T i- Refs. 7 and 8 both observed, 
through their simulation studies, that for large times t 
the correlation function G(x,t) differs from zero for dis- 
tances not too much larger than a dynamic correlation 
length £(t, T). Simulations of Marinari et al. ||7| ob- 



magnetic field through the glass temperature T g to the 
measuring temperature T, waiting a time t w , then cut- 
ting the field to zero and measuring the decay of the 
magnetization. This generates the response function, 



S(t) = d 



-M 



THAI 



H 



I dint 



(2) 



where Mtrm(1,1w) is the thermoremanent magnetiza- 
tion at time t after cutting the magnetic field to zero. 

Our approach, justified previously through magnetic 
field cycling |J and used to determine the Parisi phys- 
ical order parameter P(q), [0 makes use of the scaling 
relationship introduced by Vincent et al. Q . They show 
that barrier heights surmounted during aging are reduced 
upon a change in magnetic field by E z , a quantity related 
to the change in Zeeman energy. Our model jjj assumes 
that barriers A < E z in the initial state (before the 
magnetic field is changed) are quenched, with the popu- 
lation of occupied states transitioning essentially instan- 
taneously to those states of lowest energy corresponding 
to the new value of the magnetic field. The associated 
change in magnetization is referred to as the reversible 
part of the magnetization. The states "left behind" com- 
prise the irreversible component, and decay by diffusion 
to the sink created by the quenching of barriers A < E z . 

At small magnetic field changes, states at the barrier 
height corresponding to the waiting time t w , 



A(t w , T) = k B T(lnt w - £nr ) 



(3) 



tain satisfactory fits for £(t, T) oc (t/r ) QT / T9 , appropri- where 1/tq is an attempt frequency w ksTg/h, dif- 
ate to power law dynamics, |Q while Kisker et al. ||] fit 
satisfactorily both this proportionality and equally well 
£(t, T) oc [(T/T g )ln(t/T )] 1 /' 4 ' , appropriate to activated 
dynamics. f>| 

Our measurements consist of cooling a sample in a 



fuse towards the sink, resulting in a local equilibration 
of state occupancies on a time scale of the order of t w . 
This causes [[l0| a peak in the response function S(i), 
Eq. (2), for measurement times t close to the waiting 
time t w . When the magnetic field change is increased, 
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FIG. 1. (a)Plots of the data points for —Mtrm/H vs mea- 
surement time t, with the solid lines the analytic fit; and 
(h)S(t) defined in Eq. (2) vs t for Cu : Mn Qat.% at various 
changes in magnetic field. S(t) was calculated by taking the 
time derivative of the analytic fit (solid line) in (a). 



E z increases, and the effective barrier height from which 
diffusion takes place to the sink is A(t w , T) — E z . This 
results in a reduction of the measuring time at which S(t) 
peaks, and was first noted for experiments upon the in- 
sulating thiospinel CdCr1.-rIno.3S4 by Vincent et al, ||] 
and the amorphous system {Fe x Nii— x )j5pieBeAl3 by 
C. Djurberg et al. Jll| We follow their analysis and asso- 
ciate the peak in S(t) with an effective waiting time,i^, 
so that as a function of magnetic field change, 



A(t w , T)-E z = k B T(£ntf - tnro) 



(4) 



Systematic experiments were carried out by us for the 
metallic spin glass Cu : Mn 6at.% and for the thiospinel. 
[ p"2[ Representative curves for —Mtrm/H, S(t) are ex- 
hibited in Fig. 1. 

Studies were performed for a variety of waiting times, 
t w , magnetic field changes, H, and temperatures, T, en- 
abling us to extract values of E z over a wide range of 
parameter space. 

We interpret E z to be the magnetic energy associated 
with a change in magnetic field. In the field cooled state, 
we associate each spin to possess a magnetic susceptibil- 
ity Xfc per spin, which we can calculate from the total 
value of the field cooled magnetic susceptibility Mf c /H 
by dividing by the total number of spins N. Then the 
number of spins participating in barrier quenching (and 
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FIG. 2. A plot of log t e J s (equivalently, E. 
vs H 2 for Cu : Mn 6at.% (T/T g = 0.83, t w = 



CdCn.7lnO.3S4 {T/Tg 



from Eq. 4) 
480 sec) and 
0.72, t w = 3,410 sec) at fixed t w 
and T. The dependence is linear for magnetic fields less than 
250G and 45G, respectively, then breaks away to a slower 
dependence. The ratio of these "break" fields is in agreement 
with the respective de Almeida- Thouless lines. We do not 
have a satisfactory explanation for this change in slope. A 
different description, Ref. 2, predicts a linear dependence of 
E z upon H, which can be made to fit the data over the entire 
range of H for the thiospinel, but with a significant deviation 
at small field changes. We have chosen to focus in the text on 
the fit at small field changes because we want to be certain 
to be in the linear regime for our analysis. 



barrier hopping) can be derived from setting, 
E z = N sXfc H 2 ■ 



(5) 



The quantity A^ s defines a volume over which the spins 
are effectively locked together for barrier hopping, PJl^j 
the radius of which we define as the spin glass correlation 
length £. Were N s to be the sum of smaller independent 
clusters, the activation energy would relate to smaller 
barriers, and S(t) would not shift its peak from the vicin- 
ity of measurement times ~ t w . Extensive studies were 
made of both Cu : Mn Qat.% and CdCr1.jIno.3S4 to 
examine the H 2 dependence predicted in Eq. (5). Ex- 
perimental results for both systems, with H scales ad- 
justed to their respective de Almeida-Thouless lines, H] 
are exhibited in Fig. 2. 

Experimentally, the measured quantity N s depends 
both upon waiting time t w and temperature T. It rep- 
resents the number of correlated spin units which flip 
together within a volume £ 3 growing as a function of t w 
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FIG. 3. A plot of Ns, the number of spins participati 
in barrier quenching (and barrier hopping) from Eq. (4) i 
Cu : Mn 6at.%, as a function of Log t w at fixed temperatr 
T = 0.89T 9 = 28K. The solid curve drawn through the poii 
is the prediction for power law dynamics, Ref. 4, while t 
dashed curve is the prediction for activated dynamics, Ref. _ , 
with their exchange factor set equal to T g (i.e. independent 
of T and t). As can be seen from the two curves, the two fits 
are equally good. 



and T. 

We summarize our results for iV s in Cu : Mn 6at.% 
as a function of waiting time t w for fixed temperature 
T = 0.89T g = 28K in Fig. 3. The solid curve in Fig. 3 
is a fit to the form derived from power law dynamics, 



£(t w , T) = 0.653(Wt ) 



0.169T/T g 



(6) 



with 1/to = 4.1 x 10 12 sec 1 . A similar analysis for 
the thiospincl at T = 0.72T g = 12K, but for only two 
waiting times, 1,100 and 3,410 sec, gives £(t w , T) = 



0.53(Wr 



0.132T/T g 



The two expressions are remark- 



ably close, and suggest that Eq. (6) may be universal 
(see Fig. 4 below). The dashed curve in Fig. 3 is a fit to 
the predictions of activated dynamics, || leading to, 

at w , T) = 10- 5 [(T/T g )£n(t w /r )} 1/0 - 21 , (7) 

where we have taken the nominally t, T dependent ex- 
change factor of Ref. 5 to be a constant equal to T g . 
As seen from Fig. 3, both Eqs. (6) and (7) fit remark- 
ably well. However, the factors in Eq. (7) deserve com- 
ment. The prefactor is very troublesome, for it would 
give £ « 0(1) for t w /r a = 10 5 , while simulations @,|Jl5| 
clearly see an increasing dynamic correlation length in 
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FIG. 4. A plot of N s on a log scale vs (T/T g )£n{t w /T ) for 
Cu : Mn 6at.% (solid circles and squares) and the thiospinel 
(open triangles). The solid circles are for variable T and fixed 
t w = 8 mm, and the solid squares for variable t w and fixed 
T = 0.89T 9 = 28K for Cu : Mn 6At.%. The open trian- 



gles are for variable t v 
thiospinel. 



and fixed T = 0.72T o 



12K for the 



the regime t/ro < 10 5 . The exponent in Eq. (7), equal 
to 1/tp, is also at the lower allowable value, ||5j ip = 0.2. 

A comparison of Eq. (6) with the numerical results of 
the simulations in Refs. 0, ||, and jl5| is remarkable. 
The prefactor in the first is of order unity, with an ex- 
ponent equal to 0.16T/T g (compare with our exponent 
0.169T/T 9 ), very close to our fitted results. The second 
finds a prefactor of order unity and an exponent between 
0.12T/T g and 0.13T/T ff , and the third a prefactor close 
to unity, with an exponent 0.13T/T g . All three are in 
substantial numerical agreement with the experimental 
result, Eq. (6). 

In addition, Kisker et al. || find a proportionality of 
the same form as Eq. (7) for activated dynamics, || 
but with an exponent 1/ip = 1/0.71 ± 0.02, somewhat 
different from our fitted value 1/tp — 1/0.21. However, 
we have set the exchange factor of Ref. 5 equal to a 
constant, T g . This quantity may, however, possess a t w 
and T dependence over the temperature region of our 
measurements. Jl6|] 

In order to test the universality of our fit to the power 
law form for N s as given by Eq. (6), we plot in Fig. 
4, N s on a log scale against (T/T g )\n(t w /T ), for both 
Cu : Mn 6at.% and the thiospinel. It is seen that all of 
our data falls roughly along a universal line. It should 
be noted that two data points for the thiospinel, taken 
at low temperatures [0.36 and 0.48 xT/T g , respectively], 
would not fall on the data line in Fig. 4, and are not 



3 



exhibited. Both points correspond to data taken by a 
different method from those exhibited in Fig. 4, in that 
the peak in S(t) could not be observed, and an estimate 
of the effective age of the curves was obtained by a scaling 
procedure (see Fig. 5ofRef. 2). The precise fitting could 
be in error, but we cannot rule out that the "line" in Fig. 
4 may crossover to a weaker slope at low temperatures. 

The values of N s extracted by us from Eq. (4) are 
in the range of 10 5 , varying with t w and T according to 
Eqs. (6) and (7). This suggests that there is a meso- 
scopic number of spins participating as a unit in the ag- 
ing process, as argued initially in Ref. 2. Lederman et 
al. and Hammann et al. [|~7| suggest a large highly corre- 
lated spin system is involved in the dynamical processes 
of aging and irreversibility. 

The nature of G(x,t) in Eq. (1) lends itself to such 
an interpretation. The spin realizations Oi and Ti have 
zero overlap. Yet their correlation is finite so long as 
x < £(t, T). These realizations can be interpreted as 
rotations one from another, consistent with a change in 
Zeeman energy, relating the theoretical quantity G{x, t) 
to N s . 

The power law dynamics form for T), exhibited 

in Eq. (6), and the maximum occupied barrier A(t w , T), 
stipulated by Eq. (3), generate an interesting relation- 
ship. Substituting t w from Eq. (6) into Eq. (3) leads 
to, 
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A(t w , T)/k B T g = 6[£n^t w ,T) + 0.44] 



(8) 



This result, supported by our experimental data, was an- 
ticipated by Koper and Hilhorst Q as well as by Rieger 
jl3| and Sibani et al. Il8|. The form of Eq. (8) was also 
exhibited qualitatively in an analysis similar to our own 



by P. Granberg et al. on Cu : Mn I0at.%. JT9J Activated 
dynamics would yield a power law relationship between 
A(t w , T) and£(t w , T). 

Equation (8) shows that it would be wrong for a finite 
time scale to think of a rigid set of barriers extending to 
infinite height (the so-called "pure states"). Barriers are 
created as the correlation length increases, with detailed 
balance obeyed throughout the hierarchy. 

In summary, we have interpreted the magnetic field 
dependence of the response function to generate the 
number of spins, N s {t w , T), locked together in the 
barrier quenching (hopping) process. We identify 
[N s (t w , T)] 1 / 3 w £(t w , T) with the spin glass correla- 
tion length. The fitted expressions are nearly identical for 
measurements on Cu : Mn 6at.%, and for the thiospinel. 
This suggests that power law dynamics 0] may well be 
universal. Finally, our fitted expressions for t;(t Wl T) 
were very close to the numerical simulations of Refs. 7, 
8, and 15. 

The authors thank Drs. J. -Ph. Bouchaud, M. Mezard, 
H. Rieger, and P. Sibani for very helpful conversations 
and correspondence. We are especially indebted to Pro- 
fessor G. Parisi for pointing out the t w and T dependence 
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